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ABSTRACT 

We present an approach that improves the search for rehable astrophysical parame- 
ters (e.g. age, mass, and distance) of differentially-reddened, pre-main sequence-rich 
star clusters. It involves simulating conditions related to the early-cluster phases, in 
particular the differential and foreground reddenings, and internal age spread. Given 
the loose constraints imposed by these factors, the derivation of parameters based 
only on photometry may be uncertain, especially for the poorly-populated clusters. 
We consider a wide range of cluster (i) mass and (ii) age, and different values of (Hi) 
distance modulus, (iv) differential and (v) foreground reddenings. Photometric errors 
and their relation with magnitude are also taken into account. We also investigate 
how the presence of unresolved binaries affect the derived parameters. For each set of 
(i) - (v) we build the corresponding model Hess diagram, and compute the root mean 
squared residual with respect to the observed Hess diagram. The parameters that pro- 
duce the minimum residuals between model and observed Hess diagrams are searched 
by exploring the full parameter space of (i) - (v) hy means of brute force, which may 
be time consuming but efficient. Control tests show that an adequate convergence is 
achieved allowing for solutions with residuals 10% higher than the absolute minimum. 
Compared to a colour-magnitude diagram containing only single stars, the presence 
of 100% of unresolved binaries has little effect on cluster age, foreground and differen- 
tial reddenings; significant differences show up in the cluster mass and distance from 
the Sun. Our approach shows to be successful in minimising the subjectiveness when 
deriving fundamental parameters of young star clusters. 

Key words: ( Galaxy:) open clusters and associations: general 



1 INTRODUCTION 

Because of important dynamically-induced structural 
changes, the first few 10^ yr represent the most critical pe- 
riod in a star cluster's life, especially for the low-mass em- 
bedded clusters (ECs). At this stage, cluster dissolution is es- 
sentially related to the impulsive parental gas removal by su- 
pernovae and massive star winds. Following the rapid change 
in the gravitational potential - and the reduced escape ve- 
locity - an important fraction of the stars, even all stars, 
escape to the field. This process has been shown capable of 
dissolving most of the very young star clusters on a time- 
scale of 10 - 40 Myr (e.g. lTutukovlll978l : iGoodwin fc BastianI 
l2006f ). We point out that, because of gas and dust ejec- 
tion, young clusters with most of the stars still in the 
PMS are not n ecessarily ECs. E xamples might be B ochum 1 
llBica. Bonatto & Dutra [20081') . PismisS ( Bo natto fc Bical 
l2009bl ). and NGC 4755 lBonatto et al.ll2006l ). 

Current estimates (e.g. iLada fc Ladal l2003l : 



iBonatto fc Bical l2011bl ) suggest that less than ~ 5% 
of the Galactic ECs dynamically evolve into gravitationally 
bound open clusters (OCs). Thus, such a massive early 
dissolution of ECs may be an important source of Galac tic 
field stars (e.g. iMassev. Johnson fc Gioia-Eastwoodlll995^ . 

It is in this context that a robust determination of fun- 
damental parameters (e.g. age, distance, mass, reddening, 
etc) of young star clusters is important. However, several fac- 
tors severely challenge this task, most of which are related 
to the conspicuous presence of pre-main sequence (PMS) 
stars in Colour-Magnitude Diagrams (CMDs) of young, low- 
mass star clusters. Evolving towards the under-populated 
MS, the PMS stars are usually shrouded within a non- 
uniform dust distribution, which may lead to a high de- 
gree of differential reddening. Examples of clusters charac- 
terised by such CMDs are NGC 6611, NGC 4755, NGC 2244, 
Bochum 1, Pismi s 5, NGC 1931, vdB80, Cr 197 a nd vdB 92 
IIBonatto fc Bical l2009al : iBonatto fc Bical l2010al and refer- 
ences therein) . This effect is particularly critical for the very 
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young clusters that are still in the (gas and dust) embedded 
phase. As an additional complicating factor, star formation 
within a cluster is not characterised by a single event. In- 
stead, stars in young clus ters are observed to form over a sig- 
nificant time spread (e.g. IStauffer et al.|[l997l . and references 
therein), usually comparable to the cluster age. In summary, 
CMDs of clusters younger than tage ~ 30 — 40 Myr are ex- 
pected to contain stars with a range of ages (sC tagu) and 
affected by varying degrees of differential reddening. And, 
since most of the stellar mass of a young cluster is stored in 
the PMS, the combined effect of the age spread and differen- 
tial reddening complicates the straightforward derivation of 
cluster fundamental parameters. This is especially true when 
photometry is the only available information on a cluster. 

Given the relevance of the above issues, it's natural that 
previous approaches with different sophistication degrees 
have been develope d, exploring similar lines as in the present 
paper. For instance. iNavlor fc Jeffrie j (|2006l ) present a pow- 
erful and formally elegant maximum-likelihood method to 
derive distances and ages of young clusters through com- 
prehensive Hess diagram simulations. Although allowing for 
the presence of binaries in an isochronal population, their 
method does not include age spread or differential redden- 
ing. Also, it appears to apply more consistently to clusters 
with CMDs satisfactorily described by a single isochrone, 
i.e., those without a signifi cant PMS age spread (usually 
older than ~ 30 Myr). Later^^ da Rio. Gouliermis &: Gennard 
l|2010t ) improves on the iNavlor fc Jeffriej (|2006l ') method 
by including differential reddening, age spreads and PMS 
stars in the simulations. H owever, distance and reddening 
are no t free parameters in Ida Rio. Gouliermis fc Gennard 
instead, they are adopted from previously estimated 
values. iHillenbrand. Bauermeister fc White! (|2008l ) presents 
another attempt to modelling CMDs of young star clusters 
by means of varying star formation histories. Based on con- 
fusion between signal and noise in CMDs, they conclude that 
there is only marginal evidence for moderate age spreads in 
recent s tar forming regions an d young open clusters. More 
recentlv. lStead fc Hoard l|201ll l apply a Monte Carlo method 
to the age determination of embedded clusters through near- 
infrared (UKIDSS) photometry. They deredden the photom- 
etry of a real cluster and compare it with models built from 
theoretical isochrones. 

In this work we simulate some relevant conditions that 
usually apply to the early cluster phases to approach the 
problem of obtaining reliable fundamental parameters of 
PMS-rich clusters affected by varying degrees of differen- 
tial reddening. We adopt a simplistic approach, keeping the 
number of assumptions - and free parameters - to a mini- 
mum. In short, we start by building a distribution of arti- 
ficial stars (i.e. with mass and absolute luminosity) corre- 
sponding to clusters of a range of masses and ages. Next, 
we apply a set of values of foreground reddening, distance 
modulus and differential reddening to the model stars, and 
build the corresponding Hess diagram. Photometric uncer- 
tainties and their smearing effect on CMDs are explicitly 
taken into account. At each step we compare the artificial 
and observed Hess diagrams, searching for the set of values 
that produce the best match. Binaries - in varying fractions 
- are also included in the simulated Hess diagrams. Formally, 
our brute-force approach is not as elegant as some of the pre- 
vious methods (see above). However, it has the advantage 
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Figure 1. Top: Binary mass-ratio distribution that result from 
pairing stars with similar ages (and any mass) on a cluster un- 
dergoing continuous star formation for 20 Myr. ms and rup are 
the secondary and primary-star mass, respectively. The distribu- 
tion's peak at g fs; 0.25 has been normalized to unity. Bottom: 
J-band photometric completeness as a function of different error 
tolerances (see Sect. [2T2l for details). 



of fully exploring the parameter space in the search for the 
best set of values, although at the cost of heavy computer 
time when working with fine grids (Sect.[2]2}. 

The present paper is organised as follows. In Sect. [2] 
we discuss the relevant effects that affect the CMDs and 
describe the approach. In Sect. |3] we apply the approach 
to some control cases and discuss the results in terms of 
the residual statistics. In Sect. |4] we do the same with two 
previously studied young clusters. Concluding remarks are 
given in Sect. [5l 



2 YOUNG CLUSTER SIMULATIONS 

For consiste ncy with previous work by our group on young 
clusters (e.g. iBonatto fc Bicall20l"0al ). we consider here pho- 
tometric properties (e.g. the relationship between errors 
with apparent magnitude, quality control, CMDs) that usu- 
ally apply to 2MAS^|. The all-sky coverage and unifor- 
mity of 2MASS, associated with a moderate near-infrared 
photometric depth, provide an adequate environment for 
probing properties even of deeply embedded clusters (e.g. 



^ The Two Micron All Sky Survey, All Sky data release 
llSkrutskie et al.ll2006l) 
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Figure 2. Template CMD of a 20Myr and 250 M0 cluster with 
stars submitted to A^^ = (left panels) and A^^ = 4 (right), 
displaced at (m — M)j = (top) and (m — M)j = 11 (bot- 
tom). The arrow shows the reddening vector for Ay = 5. The 
foreground reddening is zero in all cases. For comparison, evolu- 
tionary sequences corresponding to binary fractions of 0% (filled 
symbols) and 100% (empty) are considered. 



iBonatto fc Bicall2011al and references therein). For similar 
reasons, we work with J x ( J — Ks) CMDs (this colour is 
the least affected by photo metric errors and the bes t discrim- 
inant for PMS stars - e.g. iBonatto fc Bical l2010bl ) and the 
corresponding Hess diagrams. Thus, in what follows we refer 
to the apparent distance modulus as (m — foreground 
reddening as E{J — Ks), and differential reddening as Ay^. 
Reddening transformations are based on the absorption rela- 
tions Aj/Ay = 0.276 and AK jAy = 0.118, with Ry = 3.1 
llCardelli. Clayton fc Mathid Il989l : iDutra. Santiago fc Bical 
|2002| ). 

Photometric uncertainties - which are assumed to be 
Normally distributed - are taken into account when building 
the Hess diagrams. Formally, if the magnitude (or colour) of 
a given star is given hy x ± ax, the probability of finding it 

at a specific value x is given by P{x) = ^^^^ e J , 

Thus, for each star we compute the fraction of the magnitude 
and colour that occurs in a given bin of a Hess diagram, 
which corresponds to the difference of the error functions 
computed at the bin borders. By definition, summing the 
colour and magnitude density over all Hess bins results in 
the number of input stars. As a compromise between CMD 
resolution and computational time, the Hess diagrams used 



in this work consist of magnitude and colour bins of size 
AJ = 0.2 and A(J - Ks) = 0.02, respectively. 

Binary systems are expected to survive the early evo- 
lutionary phase of low-mass clusters, with the unresolved 
pairs being somewhat brighter than the single stars and pro- 
ducing some broadenin g of the CMD stellar sequences (e.g. 
iNavlor fc Jeffrie^ |2006| ) . In addition, binaries also produce 
changes on t he initial mass functio n of young, massive star 
clusters (e.g. lWeidner. Kroupa fc Maschberger. 2009, 1. Thus, 
we also include them in our simulations by means of the 
parameter fbin, which measures the fraction of unresolved 
binaries in a CMD. According to this definition, a CMD 
having Ncmd detections, but characterised by the binary 
fraction fun, would have a number of individual stars ex- 
pressed as Af, = (1 -(- fbin)NcMD ■ For consistency with our 
assumption of continuous star formation (Sect. [22]), binaries 
are formed by pairing stars with the closest ages, regardless 
of the individual masses. This gives rise to a secondary to 
primary star mass-ratio {q — ms/mp) that smoothly in- 
creases from very-low values up to g ~ 0.25, and decreases 
for higher values of q (Fig. [l}. 



2.1 A first look at the problem 

Most of the difficulties associated with obtaining reli- 
able fundamental parameters for young clusters are sum- 
marised in Fig. [2] in which we present a template 
CMD corresponding to a 250 M0 cluster with an age 
of 20 Myr, and no foreground reddening. For the stel- 
lar mass/luminosity relati on we use the so lar-metallicity 
isoc hrone sets of Padova llCirardi et al.l 1200 2) isochronei|f| 
and ISiess. Dufour fc Forestini ( 2000l ). Both isochrone sets 
have been merged for each age considered, since Padova 
isochrones should be used only for the MS (or more evolved 
sequences), while those of Siess apply to the PMS. The merg- 
ing point occurs at the MS, at 6.5 Mq for the isochrones 
younger than 8 Myr, 5.5 M© for 10 Myr, 4.5 M© for 20 Myr, 
and 3.5 Mq for 30 Myr. The models are built with stars with 
mass ^ 0.1 Mq, the lowest mass considered in the PMS 
isochrones. 

First we consider zero distance modulus and no differ- 
ential reddening. At this point the template stars are rel- 
atively bright and simply distribute among the isochrones 
according to the age and mass, with very small photomet- 
ric uncertainties. When the template stars are displaced 
by e.g. (m — M),j — 11 (still with no differential redden- 
ing), a significant scatter shows up, increasing for fainter 
stars. Finally, when a moderate value of differential redden- 
ing {Ay^ = 4) is added to the template, the scatter tends 
to mask any relationship between stars and isochrones, even 
for (m — M)j = 0. Clearly, this effect increases with dis- 
tance modulus. At each step we also show the stellar se- 
quences corresponding to the maximum possible binary frac- 
tion, fbin = 1.0. While the mild brightening implied by the 
binaries with respect to the single stars (Sect.|4]) is clearly 
seen in all cases, the broadening ends up drowned both by 
the intrinsic PMS age spread and photometric errors. Thus, 



^ Computed for the 2MASS 
I http://stev. oapd. inaf. it/ cgi- bin/cmd\ 
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the binary-related broadening should be more noticeable in 
the MS of massive clusters of any age. 

Note that, as a photometric quality control, we have 
kept only the stars with J and Ks errors lower than 0.5. 
Such a loose constraint is used here only for illustrative 
purposes. In what follows we work only stars with errors 
lower than 0.2. The error restriction is used to somehow 
emulate the photometric completeness function of the ob- 
servations. We illustrate this effect on the luminosity func- 
tion of an artificial cluster containing an arbitrarily large 
number of stars, for statistical purposes (Fig. [T]). After fol- 
lowing the steps described below for assigning magnitudes 
and uncertainties, we applied the restriction of keeping only 
stars with aj^Q.l and 0.2. Compared to the complete lu- 
minosity function, which presents a turnover at J ~ 16.9, 
the restricted functions have turnovers that smoothly shift 
towards brighter magnitudes. 

2.2 Description of the approach 

Briefly put, we start by building the Ifess diagram corre- 
sponding to the J X (J — Ks) CMD of a young cluster. Then 
we search - by means of brute force - for the set of values of 
cluster stellar mass (Mdu), age (tage), differential redden- 
ing (Ay^), foreground reddening {E(J — Ks)), and appar- 
ent distance modulus ((m — M)j), that produces the best 
match between the observed and simulated Hess diagrams. 
Note that, to minimise the number of free parameters, we 
assume a uniform (or flat) distribution for the differential 
reddening. Alternative shapes might be tested, such as a nor- 
mal distribution around a mean value. However, this would 
require an additional parameter (the standard deviation), 
and CMDs of young clusters might lack constraints to find 
the best values for a large number of free parameters. Also, 
we assume that for a cluster of age tage and mass Mdu, stars 
(of mass m) form continuously between ^5 f 55 tage- In this 
context, the cluster age also characterises the star-formation 
time spread, since we assume the cluster age to coincide with 
the beginning of the star formation. As a caveat we note that 
assuming a flat age distribution may be somewhat unreal- 
istic, especially for clusters older than ~ 10 Myr, since this 
would imply a very slow and steady star formation rate. 
However, except for the artificial clusters used as control 
tests below, in this paper we apply our approach to clusters 
younger than this threshold. 

The approach takes on the following steps: (i) Start 



with an artificial cluster of mass Ala 



age ta 



appar- 



ent distance modulus (m — M)j, and foreground reddening 
E{J — Ks). (^ii j Sele c t the stellar masses by randomly taking 
values from lKroupal (|200ll l mass function, until the individ- 
ual mass sum yields Mdu- (Hi) Randomly assign each star 
an age ^ tage- (iv) As another simplifying assumption, con- 
sider that each star can be randomly absorbed by any value 
of (differential) reddening in the range 
(v) Apply shifts in colour and magnitude according to the 
values of E{J — Ks) and (m — M)j. (vi) Assign each arti- 
ficial star a photometric uncertainty based on the average 
2MASS errors and magnitude relationship, (vii) For more 
realistic representativeness, add some photometric noise to 
the stars. This step is done to minimise the probability of 
stars with the same mass having exactly the same observed 
magnitude, colour, and uncertainty in the CMD. Consider 



a star (of mass m and age tage) with an intrinsic (i.e., mea- 
sured from the corresponding isochrone) magnitude magi 
and assigned uncertainty <Jmag- The noise-added magni- 
tude mag is then randomly computed from a Normal dis- 
tribution with a mean magi and standard deviation amag- 
(viii) Apply the same detection limit to the model CMD 
as for the observations, so that model and data share a 
similar photometric completeness function; in practice, this 
means that stars with photometric errors higher than 0.2 
are discarded, (ix) Build the corresponding J x {J — Ks) 
Hess diagram, compare it with the observed one and com- 
pute the root mean squared residual Rrms (see below) for 
this set of values, (x) Repeat steps (i) - (x) for a range 
of values of Mdu, tage, A^^ , E{J -Ks), and {m-M)j- 
Finally, analyse the topology of the hyperspace defined 



by Rrms = Rrms[Mc, 



A'- 



,E[.J -Ks),{m-M)j) 



and search for solutions around the minimum values of 
Rr ms- Som e techn ical details are described below. 

iKroupal (|200ll ) mass function is defined as dN /dm oc 
m"<^+'^', with the slopes x = 0.3 for 0.08 m(Mo) ^ 0.5 
and X ~ 1-3 for m{MQ) > 0.5. The relation of errors 
with magnitude for 2MASS is well represented by aj = 
0.0214 + 2.48 X 10"** exp (J/1.071) and ctk, = 0.0193 -f 
9.59 X 10~® exp (/S's/1.067). As for the photometric noise 
for a star with magnitude mag ± Omag, a new magnitude is 
randomly taken from the Normal distribution characterised 
by the mean mag and standard deviation Omag- 

With respect to the root mean squared residual between 
observed {Hobs) and simulated {Hsim) Hess diagrams com- 
posed of ric and n-m. colour and magnitude bins, we define 



Rrms = 



\ 



\Hohs{i,j) - Hsim{i,j)f 



To minimise the critical stochasticity associated with 
low-mass clusters, in step (t) we build Nsim clusters of mass 
Meiu and age tage- To improve the statistical significance, 
the final simulations, in most cases, are run with Nsim = 
250. In this sense, the artificial Hess diagram (step (ix)) 
corresponds to the average density over the Nsim clusters. 

For practical reasons, the age grid is restricted to tage = 
0.2, 1, 3, 5, 8, 10, 20 and 30 Myr. Indeed, as can be in- 
ferred from Fig.O a finer age grid would be redundant. The 
stellar mass/luminosity relation is taken from the respective 
merged set of solar- met allicity Padova and Siess isochrones. 
Magnitudes for stars with intermediate age values are ob- 
tained by interpolation among the neighbouring isochrones. 
According to the adopted isochrone sets, the minimum stel- 
lar mass is 0.1 M©, while the maximum ranges from 60 M© 
(at 0.2Myr), 36M0 (5Myr), 19Mq (lOMyr), and 9M0 
(30 Myr). 

As described above, the brute-force nature of our ap- 
proach tends to be very time consuming. For instance, a 
typical simulation for a cluster with ~ 500 stars would re- 
quire a parameter grid composed (at least) of 21 mass bins, 
5 ages, 21 distance moduli, 21 foreground reddenings and 21 
differential reddenings. Then, including the Nsim clusters in 
the simulation, the runtime (for the minimum grid) is about 
10 hours on a single core of an Intel Core i7 920@2-67 GHz 
processor. 
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Figure 3. Rrms contour maps for selected parameters of 
MODEL#2 (left) and MODEL#3 (right). Darker colours indicate 
lower values of Rrms- The absolute minima are Rrms = 0.034 
(MODEL#2) and Rrms = 0-028 (MODEL#3). 



Figure 4. Observed (top panels) and simulated (bottom) Hess 
diagrams for MODEL#2 (left) and MODEL#3 (right). Lighter 
colours indicate higher densities of stars. The simulated diagrams 
have been built with the 10%-higher parameters (Table [TJ. 



3 CONTROL EXPERIMENTS 

Before turning to actual cases, we apply the approach de- 
scribed above to template CMDs built with typical param- 
eters found in PMS-rich young clusters. The main reason 
is to examine the ability to recover the input parameters 
and their relation with the Rrms statistics. The relevant 
parameters of the adopted models cover the ranges 100 ^ 
M^iuiUe) s; 500, 3 < tageiMyr) 20, 9.5 ^ {m-M)j ^ 
12, 0.1 ^ E{J - Ks) sS 0.4, and 0.5 ^ ^ 3.5. The mod- 
els are described in Table [T] which also contains quantitative 
details of the search for solutions on the Rrms maps. 

Besides the single solution for the absolute Rrms min- 
imum, we also explore the solutions with Rrms occurring 
within thresholds 5%, 10%, 25%, and 50% higher than the 
absolute minimum. For these, the parameters given in Ta- 
ble [T] correspond to the weighted average of all the solutions 
matching each threshold. As weight for each solution we 
take the individual value of 1/R?rms- Except for a few cases, 
there's little difference among the average values of a given 
model, from the absolute minimum to the 50% threshold. 
As expected, the most noticeable difference in the output 
lies in the significantly increasing dispersion around the av- 
erage for higher thresholds. Although somewhat subjective, 
we believe that the best convergence is obtained with the 10% 



threshold (which also produces realistic errors) , which is con- 
sistent with the somewhat irregular topology of the Rrms 
maps (see below). Overall, the approach seems very sensi- 
tive to the parameters, especially the age and differential 
reddening. The residual Rrms of the 10%-higher solutions 
occur in the range 0.023 Rrms ^ 0.051. As an additional 
perspective on the statistical significance of the solutions 
within the adopted residual thresholds, we also provide in 
Table[l]the percentage of Nsoi with respect to the full range 
of possibilities (which depends on the number of parameter 
bins). Note that even allowing for the 10%-higher threshold, 
the fraction of acceptable solutions corresponds to ~ 1 — 5% 
of the total number. 

Selected two-dimensional projections of Rrms ~ 
Rrms{McIu, tage, Ay^ , E{J - K s) , (m - M) j) are shown in 
Fig. 13] for M0DELS#2 and 3, those that present extreme 
values of cluster mass and differential reddening. As dis- 
cussed above, the approach presents conspicuous conver- 
gence towards the input values of the age, cluster mass, 
and apparent distance modulus. However, the Rrms maps 
present some spread in both foreground and differential 
reddening values. In particular, there clearly is an anti- 
correlation between both reddening sources, in the sense 
that low (or high) foreground reddenings are compensated 
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Table 1. Recovery of model fundamental parameters (Ngi^ = 500; fun = 0.0) 



Rrms 
(1) 


Range 
(2) 


(3) 


fsol 

{%) 
(4) 


Mctu 

(Mq) 
(5) 


Age 
(Myr) 
(6) 


(m - M)j 
(mag) 
(7) 


E{.J-Ks) 
(mag) 
(8) 


aDR 
(mag) 
(9) 


McMD 
(Mq) 
(10) 




Input parameters 


of MODEL#l 


250 


3 


9.5 


0.1 


3.0 




0.0463 


Abs. minimum 


1 


1.7 X 10-3 


240 ± 


3.0 ±0.0 


9.60 ±0.00 


0.14 ±0.00 


2.8 ±0.0 


160 ±0 


^ 0.0487 


5% higher 


303 


5.1 X 10-1 


233 ± 11 


3.0 ±0.1 


9.59 ±0.04 


0.14 ±0.02 


2.8 ±0.2 


156 ±8 


s; 0.0510 


10% higher 


1030 


1.7 


235 ± 16 


3.0 ±0.1 


9.58 ±0.05 


0.13 ±0.02 


2.9 ±0.2 


156 ± 10 


< 0.0579 


25% higher 


4002 


6.7 


234 ± 20 


3.0 ±0.1 


9.57 ±0.08 


0.13 ±0.03 


2.9 ±0.2 


155 ± 13 


< 0.0695 


50% higher 


9023 


15.0 


234 ± 21 


3.0 ±0.2 


9.56 ±0.11 


0.12 ±0.03 


2.9 ±0.2 


155 ± 14 



Input parameters of M0DEL#2: 100 5 10.0 0.2 0.5 

0.0336 Abs. minimum 1 1.4 X 10"^ 100 ± 5.0 ± 0.0 10.10 ±0.00 0.22 ± 0.00 0.2 ± 0.0 106 ± 

^ 0.0353 5% higher 214 3.0 X 10"^ 104 ± 6 5.0 ±0.1 10.08 ± 0.04 0.21 ± 0.02 0.4 ±0.2 103 ± 6 

< 0.0369 10% higher 828 1.2 104 ± 9 5.0 ±0.1 10.07 ± 0.06 0.20 ± 0.02 0.5 ± 0.3 101 ± 9 

^ 0.0420 25% higher 8977 13.0 108 ± 14 6.3 ± 1.5 9.95 ±0.17 0.20 ± 0.03 0.5 ±0.3 106 ± 14 

^ 0.0504 50% higher 28064 40.0 106 ± 15 6.5 ± 1.5 9.95 ±0.19 0.21 ± 0.03 0.5 ±0.3 103 ± 15 





[nput parameters 


of MODEL#3: 


500 




10 


12.0 


0.3 


3.5 




0.0279 


Abs. minimum 1 


5.1 X 10-"* 


495 ± 


8 


0±0.0 


12.20 ±0.00 


0.30 ±0.00 


3.7 ± 0.0 


261 ±0 


^ 0.0293 


5% higher 


1811 


9.2 X 10-1 


494 ± 22 


8 


4± 0.8 


12.16 ±0.09 


0.31 ±0.02 


3.6 ±0.2 


265 ± 11 


^ 0.0307 


10% higher 


5871 


3.0 


492 ± 23 


8 


8± 1.0 


12.12 ±0.13 


0.31 ±0.02 


3.6 ±0.2 


265 ± 12 


^ 0.0349 


25% higher 


21327 


11.0 


490 ± 24 


8 


8± 1.2 


12.11 ±0.17 


0.31 ±0.03 


3.6 ±0.2 


263 ± 13 


^ 0.0418 


50% higher 


51482 


26.0 


489 ± 24 


8 


8± 1.5 


12.12 ±0.20 


0.30 ±0.04 


3.6 ±0.2 


263 ± 13 





[nput parameters 


of MODEL#4: 




150 


20 


11.0 


0.4 


1.5 




0.0209 


Abs. minimum 1 6.1 X 10 


120 ± 


20.0 ±0.0 


11.10 ±0.00 


0.46 ± 0.00 


1.1 ± 0.0 


94 ±0 


^ 0.0219 


5% higher 


2154 


1.3 


142 ± 19 


24.3 ±4.9 


11.09 ±0.06 


0.44 ±0.03 


1.5 ±0.3 


102 ± 14 


s: 0.0230 


10% higher 


8735 


5.3 


142 ± 20 


24.5 ± 5.0 


11.09 ±0.10 


0.44 ±0.03 


1.5 ±0.3 


101 ± 14 


s; 0.0261 


25% higher 


43299 


26.0 


138 ± 23 


24.6 ±5.5 


11.07 ±0.19 


0.43 ±0.04 


1.5 ±0.3 


99± 16 


^ 0.0313 


50% higher 


89775 


55.0 


135 ± 24 


24.3 ±6.1 


11.07 ±0.22 


0.42 ±0.05 


1.5 ±0.3 


97± 17 



Table Notes. Cols. (1) and (2): Rrms value and corresponding range with respect to the absolute minimum; Col. (3); number of solutions 
occurring within the Rrms range; Col. (4): percentage of N^^i with respect to the full range of solutions; Col. (5): actual cluster mass; 
Col. (9): differential reddening; Col. (10): mass detected in the CMD. Nsim is the number of simulated clusters of mass and age 

tage- The average stellar mass of the models is m,* 0.6 Mq. 



for by the approach with high (or low) values of difTerential 
reddening. 

An important information that in principle could be ob- 
tained from CMDs is the cluster mass. More specifically, the 
fraction stored in MS and PMS stars (in this work it means 
the mass sum for all stars more massive than ^ 0.1 M©). 
However, as discussed in Sect. [U the intrinsic age spread 
together with the presence of differential reddening, compli- 
cate the task of finding the mass of each star in a CMD. 
Besides, because of limitations inherent to any photomet- 
ric system (which, in the present work, are simulated by 
means of the quality criterion (step (viii) in Sect. 12.21) . the 
number of stars that remain detectable in a CMD tends to 
decrease as the distance modulus increases. Consequently, 
the same applies to the stellar mass that would be mea- 
sured in a CMD (Mcmd) with respect to the actual cluster 
mass (Mciu)- After derivation of the fundamental parame- 
ters, Mcmd can be estimated by finding the probable mass 



for each star in the CMD. This can be done by interpolation 
(of the observed colour and magnitude of each star) among 
the nearest isochrones (for instance. Fig. [7|, but the final 
result would depend heavily on the amount of differential 
reddening, photometric noise, etc. Obviously, the presence 
of a large fraction of binaries would produce low values of 
Mcmd- By construction, the present approach provides di- 
rectly the actual cluster mass and, since we keep track of the 
mass (and photometry) of each star that is used in the sim- 
ulations, we can also compute the mass present in a CMD. 
Both mass values are given in Table [1] Interestingly, for the 
range of (m — M)j covered by the models, the actual clus- 
ter mass and that present in the CMDs are essentially the 
same. 

Another quality assessment criterion is provided by the 
compared morphology between the model and simulated 
Hess diagrams (Fig. |4|. The latter have been constructed 
with the parameters found with the 10% threshold (Table[T]). 
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Figure 5. Same as Fig. [3] for the young clusters CoIIinder 197 
(left) and PismisS (right). The absolute minima are Rums = 
0.054 (CoUinder 197) and Rrms = 0.019 (PismisS). 



Figure 6. Same as Fig. |4] for the field-star decontaminated (top 
panels) and simulated (bottom) Hess diagrams of CoUinder 197 
(left) and Pismis5 (right). 



Given the low-mass nature of the models considered here, 
some discreteness in their Hess diagrams is expected. On 
the other hand, the simulated Hess diagrams correspond to 
the average over Nsim simulated clusters (of mass Mdu and 
age tage), and thus, they present a smoother distribution. 
Nevertheless, model and simulated Hess diagrams present a 
good correspondence. 



4 APPLICATION TO ACTUAL YOUNG 
CLUSTERS 

Having demonstrated the convergence efficiency of our ap- 
proach with artificial clusters (Sect. [3]), we now move 
on to examining properties of actual cases. For this we 
have selected two young clusters previously studied by our 
group, CoUinder 197 ([B onatto & Bica 2010aj) and Pismis 5 
l|Bonatto fc Bical l2009t< ). Both present typical CMDs of 
young clusters, with the difference that CoUinder 197 has 
about 10 times more stars (essentially PMS) in the CMD 
than Pismis 5. 

To put the results derived in the present paper in con- 
text, we provide here a brief explanation of the previous 
method used by our group to obtai n parameters of both clus - 
ters. While the proc e dure u sed in iBonatto fc Bical (|2010al ) 
and iBonatto fc Bical (|2009bh to construct the field-star de- 



contaminated CMDs follows a quantitative approach, the 
fundamental parameter derivation used in both cases is 
somewhat subjective, depending essentially on a qualita- 
tive assessment of the differential reddening. Specifically, the 
MS-I-PMS isochrones (for a range of ages) are set to zero 
distance modulus and foreground reddening, and shifted in 
magnitude and colour until they produce a satisfactory fit of 
the blue border (and redwards spread) of the MS and PMS 
stellar distribution. Since we did not dispose of any quanti- 
tative tool to evaluate the differential reddening, we simply 
assumed it to be equivalent to the size of the reddening vec- 
tor that matched the colour spread among the lower-PMS 
sequences. For the (CMD extracted from t he) region within 
R = 10 ' from the centre of CoUinder 197, IBonatto fc Bical 
(|2010ah find tage = 5 ± 4Myr, (m - M)j = 10.4 ± 0.4, 
E{J-Ks) = 0.17 ±0.08, McMD ~ 450 ±100Mq, and 
the distanc e from the Sun dp, = 1. 1 ± 0.2 kpc. For R — 6' 
of Pismis 5, IBonatto fc Bical l|2009bl ) find tage = 5 ± 4 Myr, 
{m-M)j = 10.4 ± 0.1, E{J-Ks) = 0.20 ± 0.02, and 
McMD « 58 ± 8M0, and d© = 1.0 ± O.lkpc. The adopted 
age (and uncertainty) is simply the average (and half 
the difference) between the youngest and oldest MS-I-PMS 
isochrones compatible with the CMD morphology. The mass 
was estimated by multiplying the number of CMD stars by 
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Table 2. Fiindiuiu^iilal ])arainotorR of the young cliistcTs CoUindor 197 and Pismis 5 



Rrms 


Range 




fsol 






Age 


(m - M)j 


E{J - Ks) 




McMD 


d© 








(%) 




(Mo) 


(Myr) 


(mag) 


(mag) 


(mag) 


(Mo) 


(kpc) 


(1) 


(2) 


(3) 


(4) 




(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


CoUindcr 197 (fun = 0.0) 


0.0546 


Abs. min. 


1 


3.6 X 10- 


-4 


390 ±0 


10.0 ±0.0 


10.30 ± 0.00 


0.18 ±0.00 


3.1 ±0.0 


245 ±0 


1.0 ±0.0 


^ 0.0574 


5% higher 


186 


6.7 X 10" 


-2 


392 ± 18 


9.8 ±0.6 


10.37 ±0.09 


0.18 ±0.02 


3.0 ±0.2 


249 ± 11 


1.0 ±0.1 


^ 0.0601 


10% higher 


759 


2.7 X 10- 


-1 


393 ± 24 


9.5 ±0.9 


10.41 ±0.12 


0.18 ± 0.03 


3.0 ±0.3 


248 ± 15 


1.0 ±0.1 


^ 0.0683 


25% higher 


6527 


2.4 




388 ± 36 


9.8 ±2.7 


10.42 ±0.22 


0.18 ±0.04 


3.0 ±0.5 


245 ± 23 


1.0 ±0.1 


^ 0.0819 


50% higher 


28754 


10.0 




381 ± 41 


9.9 ±4.0 


10.48 ± 0.35 


0.18 ±0.05 


3.0 ±0.5 


241 ± 26 


1.1 ±0.2 


Collinderl97 (hin = 0.5) 


0.0566 


Abs. min. 


1 


2.0 X 10- 


-4 


660 ±0 


10.0 ±0.0 


10.80 ± 0.00 


0.18 ±0.00 


3.1 ±0.0 


351 ±0 


1.2 ±0.0 


^ 0.0595 


5% higher 


354 


6.9 X 10- 


-2 


666 ± 30 


10.0 ±0.2 


10.85 ± 0.07 


0.19 ±0.02 


3.1 ±0.2 


354 ± 16 


1.3 ±0.1 


^ 0.0623 


10% higher 


2037 


4.0 X 10- 


-1 


661 ± 39 


9.9 ±2.3 


10.88 ±0.21 


0.19 ± 0.03 


3.1 ±0.3 


354 ± 21 


1.3 ±0.1 


^ 0.0708 


25% higher 


22660 


4.4 




648 ± 53 


10.9 ±4.7 


10.86 ± 0.42 


0.19 ±0.04 


3.1 ±0.5 


347 ± 28 


1.3 ±0.2 


^ 0.0849 


50% higher 


84596 


17.0 




645 ± 57 


10.6 ±5.2 


10.96 ±0.49 


0.19 ±0.05 


3.0 ±0.6 


345 ± 31 


1.3 ±0.3 


Collinderl97 {fun = 1-0) 


0.0553 


Abs. min. 


1 


2.0 X 10- 


-4 


820 ±0 


10.0 ±0.0 


11.20 ±0.00 


0.20 ± 0.00 


3.0 ±0.0 


405 ±0 


1.5 ±0.0 


^ 0.0580 


5% higher 


359 


7.0 X 10- 


-2 


824 ± 35 


9.6 ±0.8 


11.21 ±0.11 


0.19 ± 0.02 


3.1 ±0.2 


404 ± 17 


1.5 ±0.1 


^ 0.0608 


10% higher 


1697 


3.3 X 10- 


-1 


816 ± 45 


9.3 ± 1.0 


11.27 ±0.14 


0.19 ±0.02 


3.1 ±0.3 


402 ± 22 


1.5 ±0.1 


^ 0.0691 


25% higher 


12174 


2.4 




803 ± 55 


9.9 ±3.1 


11.25 ±0.27 


0.19 ±0.04 


3.1 ±0.5 


395 ± 27 


1.5 ±0.2 


^ 0.0829 


50% higher 


54465 


11.0 




796 ± 58 


10.5 ±4.8 


11.27 ±0.45 


0.20 ± 0.05 


3.1 ±0.5 


390 ± 28 


1.5 ±0.3 


Pismis 5 {fun = 0.0) 


0.0188 


Abs. min. 


1 


3.5 X 10- 


-4 


60 ±0 


5.0 ±0.0 


10.75 ± 0.00 


0.08 ± 0.00 


2.9 ±0.0 


37 ±0 


1.3 ±0.0 


!g 0.0197 


5% liigher 


1301 


4.5 X 10- 


-1 


61 ±6 


5.0 ±0.1 


10.75 ± 0.08 


0.09 ± 0.03 


2.9 ±0.3 


38 ±4 


1.3 ±0.1 


^ 0.0207 


10% higher 


4816 


1.7 




62 it 10 


5.0 ±0.9 


10.79 ±0.14 


0.09 ± 0.03 


2.9 ±0.3 


39 ±6 


1.3 ±0.1 


^ 0.0235 


25% higher 


51006 


18.0 




68 ±18 


5.6 ±1.9 


10.88 ±0.26 


0.11 ±0.04 


2.9 ±0.3 


42 ±11 


1.4 ±0.2 


^ 0.0282 


50% higher 


133571 


46.0 




67 ±21 


5.8 ± 1.9 


10.90 ±0.31 


0.11 ±0.04 


2.9 ±0.3 


41 ±13 


1.4 ±0.2 


Pismis 5 {fun = 0.5) 


0.0176 


Abs. min. 


1 


3.8 X 10- 


-4 


110 ±0 


4.0 ±0.0 


11.75 ±0.00 


0.10 ±0.00 


3.0 ±0.0 


57 ±0 


2.1 ±0.0 


^ 0.0185 


5% higher 


1395 


5.3 X 10- 


-1 


111 ± 8 


4.0 ±0.2 


11.80 ±0.09 


0.10 ± 0.02 


3.0 ±0.3 


58 ±4 


2.1 ±0.1 


^ 0.0194 


10% higher 


6156 


2.3 




110 ± 13 


3.7 ±0.7 


11.89 ±0.23 


0.10 ±0.03 


3.1 ±0.3 


56 ±6 


2.2 ±0.2 


^ 0.0220 


25% higher 


52449 


20.0 




119 ± 22 


4.6 ± 1.4 


11.82 ±0.29 


0.11 ±0.04 


3.1 ±0.3 


62 ±11 


2.1 ±0.3 


^ 0.0264 


50% higher 


140800 


54.0 




121 ± 23 


5.3 ±2.0 


11.82 ±0.31 


0.11 ±0.04 


3.1 ±0.3 


58 ±11 


2.1 ±0.3 


Pismis 5 {ftin = 1.0) 


0.0189 


Abs. min. 


1 


3.2 X 10- 


-4 


140 ±0 


5.0 ±0.0 


11.60 ±0.00 


0.10 ± 0.00 


3.1 ±0.0 


67 ±0 


1.9 ±0.0 


^ 0.0198 


5% higher 


4865 


1.6 




131 ± 14 


4.3 ± 1.0 


11.74 ±0.20 


0.10 ±0.02 


3.1 ±0.3 


62 ±7 


2.0 ±0.2 


s; 0.0208 


10% higher 


20900 


6.7 




131 ± 18 


4.1 ± 1.0 


11.81 ±0.26 


0.10 ±0.03 


3.2 ±0.3 


62 ±9 


2.1 ±0.3 


^ 0.0236 


25% higher 


95284 


30.0 




134 ± 23 


4.6 ± 1.7 


11.84 ±0.32 


0.11 ±0.04 


3.2 ±0.3 


63 ±11 


2.1 ±0.3 


^ 0.0283 


50% higher 


183790 


59.0 




132 ± 24 


4.9 ±2.0 


11.88 ±0.34 


0.11 ±0.04 


3.2 ±0.3 


61 ± 11 


2.2 ±0.3 



Table Notes. CoUinder 197 presents 690 stars in the CMD, while Pismis 5 shows only 101. CoUinder 197 was analysed with Ngim = 250 
simulated clusters and, for equivalent statistical results, Pismis 5 was analysed with Ngim = 1000. is the binary fraction. The average 
stellar mass of the model clusters is mUT » 0.6 M©. 
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3 Myr 

5 Myr 

10 Myr 
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• Observed 
0.2 Myr 
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0.5 1 
(J-KJ 



1.5 



0.5 1 

(J-Ks) 



1.5 



Figure 7. Observed CMDs of Cr 197 (left) and Pismis 5 (right) 
compared to a single realisation taken from the respective simu- 
lations. The extraction radii are 10' (Crl97) and 6' (Pismis 5). 



the average stellar mass (based on lKroupalboOll MF); thus, 
it represents the CMD mass. 

Considering the above features, both clusters are excel- 
lent candidates to apply the present approach. Given the low 
number of stars (101) in the CMD of Pismis 5, we increased 
the number of simulated clusters to Nsim = 1000 for statis- 
tically more significant results. Collinder 197 has 690 stars 
in the CMD, so, we kept Nsim = 250. For a more compre- 
hensive analysis, we have applied our approach considering 
the binary fractions fbin = 0, 0.5 and 1. The fundamental 
parameters obtained with the present approach are given in 
Tabled 

Most of the parameters are essentially insensitive to the 
presence of unresolved binaries in the CMD. This is espe- 
cially true for the age, foreground reddening and, to a lesser 
degree, differential reddening. On the other hand, the most 
remarkable difference between the presence of (CMD unre- 
solved) binary fractions of 0% and 100% lies in the cluster 
mass, with fbin = 1 requiring ~ twice more mass than for 
fbin = 0. Next, the binary-related brightening of the stel- 
lar sequences (for fbin = 1) requires an apparent distance 
modulus, on average, 0.84 and 0.96 higher than for fbin ~ 0, 
respectively for Collinder 197 and Pismis 5. This reflects on 
distances to the Sun ~ 50% higher if fbin = 1- Within the 
uncertainties, the parameters corresponding to fbin = 0.5 
tend to be closer to those obtained with fbin = 1- Also, 
for a given binary fraction, the parameters are essentially 
unchanged with respect to the Rrms ranges considered, al- 



though with an increasing dispersion around the average for 
high Rrms- Finally, the convergence level of the approach, 
as measured by Rrms, is essentially the same for the bi- 
nary fractions considered here. In this sense, our approach 
appears to be insensitive to the binary fraction, at least for 
clusters with a significant age spread. 

Irrespective of the binary fraction, our results confirm 
that Pismis 5 is indeed a very-low mass {Mdu ~ 60 — 
140 M0) and young {tage ~ 4 — 6 Myr) cluster, affected by 



a moderate amount of differential reddening {A-^, 



3.0), 



and a low foreground reddening {Ay ~ 0.6). Such a low 
Av, together with the moderate apparent distance modulus 
((m — M)j ~ 10.8, for fbin = 0), puts Pismis 5 at a distance 
from the Sun of = 1.3 ± 0.1 kpc. However, this distance 
may be significantly higher (d© = 2.1 ± 0.3 kpc) if fbin = 1- 

As anticipated, Collinder 197 is more massive {Mdu ~ 
390 - 820 M0) and somewhat older {tage ~ 9.5 Myr) than 
Pismis 5, also with a similar value of differential reddening 
{Ay^' ~ 3.0), and a somewhat higher foreground reddening 
(Av ~ 1.2). The distance moduh (m - Af)j ~ 10.4 {fbin = 
0) and (m — M)j ~ 11.3 {fbin = 1) imply the distances 
cIq = 1.0 ±0.1 kpc and (Iq = 1.5 ±0.1 kpc, respectively. Un- 
der the assumptions adopted in Sect. 12.21 the star- formation 
age spread of Collinder 197 is about twice that of Pismis 5. 
Interestingly, despite the subjectiveness of the (previous) 
qualitative assessment we applied to Collinder 197 and Pis- 
mis 5, the results of both methods are somewhat compatible, 
within the uncertainties. 

The Rrms projections of Collinder 197 and Pismis 5 
(Fig. [5| present similar convergence patterns as those of the 
control tests (Fig. [31), including the foreground and differ- 
ential reddening anti-correlation. The residual Rrms of the 
10%-higher solutions of Pismis 5 {Rrms ^ 0.021) is even 
lower than those of the control tests, and somewhat higher 
{Rrms ^ 0.06) for Collinder 197. However, we point out 
that part of this may be linked to the fact that Collinder 197 
has ~ 7 times more stars in the CMD, which can result in a 
larger sum of the residuals. Next, considering the low num- 
ber of free parameters used by our approach, the observed 
and simulated Hess diagrams (Fig. [5| present a satisfactory 
correspondence. Finally, in Fig. [7|we compare a single CMD 
realisation - randomly selected among the Nsim simulated 
clusters - with the observed CMDs. Having in mind the fact 
that both are low-mass, young clusters, some differences are 
more evident because of the low number of stars, especially 
in the MS, when an isolated, random CMD is used to illus- 
trate the simulations. This effect should be particularly no- 
ticeable in poorly-populated clusters, such as Pismis 5. Nev- 
ertheless, simulated and observed CMDs are similar, in both 



5 SUMMARY AND CONCLUSIONS 

In this paper we describe an approach based on CMDs and 
near-infrared photometry to obtain more accurate funda- 
mental parameters of young star clusters. However, since the 
presence of large fractions of circumstellar disks around PMS 
stars can lead to significant excesses in the near-infrared (e.g. 
iMever. Calvet &: Hillenbrandlll997l ). the applicability of our 
appr oach may be restricted to c lusters older than ~ 3 Myr 
(e.g. iHaisch. Lada &: La dal l200ll ). Given its statistical na- 
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ture, our method should be very efficient for weU populated 
clusters, which give rise to smoother stellar-density distri- 
butions in the Hess diagrams. But more interestingly, it is 
expected to work also with the low-mass clusters with CMDs 
dominated by PMS stars and significantly affected by dif- 
ferential reddening. Among these parameters, cluster mass 
and age are important to understand the dynamical state 
(e.g. iBonatto fc Bical l2005'l of clusters, especially those un- 
dergoing the rapidly changing and potentially-destructive 
evolutionary phase characterised by the first few lO'^ yr. 

The approach involves simulating the effect of (random) 
differential reddening on MS and PMS stars, and its observ- 
able results on CMDs. In this paper we apply it to CMDs of 
star clusters simulated in the 2MASS photometric system. 
Obviously, it can be adapted to any photometric system, 
provided isochrones for low-mass stars and very young ages 
are available. 

Beginning with a given cluster mass {Mdu) and age 
(tage), wB distribute th e individual ste llar masses and magni- 
tudes according to the lKroupal (|200ll ) MF and an MS-I-PMS 
isochrone set. Subsequently, we build the corresponding sim- 
ulated Hess diagram for different values of the apparent 
distance modulus ((m — M)j), differential (Ay^) and fore- 
ground {E{J — Kb)) reddening values. We repeat this for a 
range of cluster mass and age, searching for the best match 
between the simulated and observed Hess diagrams. The 
best-match parameters are searched around the minima of 
the hypersurface defined by the root mean squared residual 

RrMS = RRMs{Mclu,tage, A^^ , E{J - K s) , {m - M) j) . 

Tests with model clusters containing parameters typi- 
cal of objects undergoing such an early phase have shown 
that the present tool converges to the input values, espe- 
cially when one allows for solutions that occur with Rrms 
10% higher than the absolute minimum. We also investigate 
how unresolved binaries affect the derived parameters. Com- 
pared to a CMD containing only single stars, even assuming 
a presence of 100% of unresolved binaries has little effect 
(to within the uncertainties) on cluster age, foreground and 
differential reddenings. Significant differences occur in the 
cluster mass and distance from the Sun. About twice more 
mass (in individual stars) is required for ftin ~ 1 than for 
fbin = 0, and because of the relative brightening due to 
the binaries, the distance from the Sun is « 50% larger for 
ftin = 1 than for /f,i„ = 0. 

As a caveat we note that we consider here a partic- 
ular isochrone set. Thus, the results of our approach tend 
to be model dependent, since different isochrone sets may 
lead to different values of mass and age of ind ividual stars 
(e.g. iHillenbrand. Bauermeister fc Whitel l200^1 and, conse- 
quently, star clusters. In addition, we minimise the number 
of free parameters by assuming uniform (or flat) distribu- 
tions of differential reddening and stellar age. Both con- 
ditions may be partly unrealistic, especially the latter for 
clusters older than ~ 10 Myr, since it would imply a slow 
and steady star formation rate. However, the clusters we 
use here to test the approach, CoUinder 197 and Pismis 5, 
are younger than ^ 10 Myr. The relative low-mass nature of 
both clusters may imply CMDs lacking constraints to find 
the best values for a large number of free parameters. 

The general conclusion is that the inclusion in the sim- 
ulations of several effects that affect the CMD morphology - 
especially the differential reddening - produces constrained 



values of the fundamental parameters. Thus, when photome- 
try is the only available information, our approach minimises 
the subjectiveness associated with the parameter derivation 
of young clusters. 
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